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Recant PtoRress Summary 


This progress report sumnarizes the recent and projected efforts in 
investigating the reduced order model problem in distributed parameter 
systems adaptive Identification and control under NASA Grant NAG**I'>7 
sponsorship. A lengthy siSHilation study 11] of the reduced order problem 
in scalar adaptive control of lumped-parameter systems, projected in a 
previous interim report, has only recently been completed. A compre- 
hensive examination (compiled over the past several months) of real-time 
centralized adaptive control options for flexible spacecraft is provided 
in the remainder of this section of this progress report. (The real-time 
objective I as used here, excludes the possibility of separating identifi- 
cation and control in time as suggested in [2],) This overview prompts 
the departure from the anticipated narrow focus on the NASA Langley beam 
control experiment and a shift to tlevolopmenL of an original, general 
approach to this problem as projected in suction II. Section III lists 
the references cited In the first two sections. Sections IV and V provide 
a listing of recent presentations and publications of work sponsored by 
NASA Grant NAG- 1-7. The final section is a cumulative list of sponsored 
papers, which have appeared in the open literature. 
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1. Introduction to Four ApproticheB to Adaptive Control of Flexible 
Spacecraft 

The fo'ir approachea to adaptive control of flexible spacecraft dis- 
cussed in the next four subsections are: 

(i) Assume that sufficiently accurate elgenshapes are provided a priori 
and perform simultaneous adaptive modal identification and control as, 
e.g., in (3) and [4]. a 

(il) Assume that eigenshapeo can be accurately approximated by a finite 
dimensional, linear combination of preselected orthogonal spatial functions 
(as suggested in (2]) and simultaneously estimate the parameters forming 
the elgenshapes and the parameters in their dynamic amplitude behavior. 
Combine tills Identification with on-line solutica of the decoupled control 
problems as In [3j and [4 ].d 

(ill) Treat the actuator/distributed parameter system (DPS)/aenaor com- 
bination as a multi-input, multi-output (MIMO) system with finite (but 
large) "state” dimension. Adaptively observe [5] [6] or identify [7] this 
MIMO system and solve the coupled MIMO control problem on-line using these 
parameter (and state) estimates in a state feedback [6, sect. 6.3] (5] or 
transfer function configuration (8, sect. 7.3] [91.Q 

(iv) Select an adequacely dimensioned feedback control structure for the 
MIMO actuator /DPS/ sensor system and directly update the controller para- 
meters to asymptotically achieve pole placement (as proposed in the scalar 
case in [10] and [11]), o 

These four approaches (especially the last three, which remain open develop- 
ment issues) are discussed in sufficient detail to pinpoint their respective 
llmltatiuns requiring further study. 
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2. Slroulcitneous Modal IdentlflGatioo arid Control 

In (3] and {4} a nodal adaptive control strategy for distributed 
parameter syetems is developed i which relies on a priori specification of 
the decoupling, spatial eigenbasls. As is commonly acknowledged, the most 
accurate modal model synthesis procedures yield eigenshape predictions 
with possibly high inaccuracy Increasing with Che spatial frequency of the 
elgenshapes. It is expected that in certain cases such a priori eigenshape 
errors in the strategy of [3] and [4] would lead to instability. The un- 
answered question is the problem dependent on: How inaccurate can these pre- 
speclfled elgenshapes be before such unacceptable behavloz results? This 
does not even consider the problems introduced by the reduced-order effects 
of modal expansion truncation as noted in (4) and [12]. The next three 
subsections are aimed at circumventing the roquireraent of exact eigenshape 
prespeclflcatlon. The study of reduced order effects would follow as noted 
in Section II. See the appendix in subsection 6 for a brief summary of 
the separated variable technique of modeling distributed parameter systems 
used in [3] and [4]. 


3. Slwulctineom Eigeniihaip* •nd Modal Dvnamica Estlwatlon 

Consider the distributed parameter system output described, as In 
{3J and {4], by the sum of products of modal elgenshape and decoupled 
amplitude dynamics 

N 

y(x^,k) - T^ti(Xq)y^Oc) (3-1) 

where y(x ,k) represents the output measured at the qth sensor location x 
at time k, ♦^(Xq) Is the magnitude of the 1th elgenshape at location x^, 
and y^(k) la the amplitude of the 1th mode at time k. Since the modal 
dynamics are uncoupled 

n 

YlOc) - I tai,iyi(k-l) ■¥ b^j^u^(k-l)], (3-2) 

I**! 

where the order of the dynamics n Is typically 2 for the linearised small 

amplitude motion of flexible structures and u^(k-S,) the 1th modal Input at 
time k-1 

c 

Ui(k) - I ♦^^(x.)u(x.,k), (3-3) 

1 j.l 1 3 J 

where u(Xj,k) la the actual force applied through the actuator located at Xj. 
Assume, as suggested In [2], that the elgenshapes can be formed as a 
flnlte-dlmenslonal, linear combination of preepeclfled Independent spatial 
functions fg(x) as 

P 

<^j^(x) - c, f (x), (3-4) 

8*1 

The simultaneous elgenshape and modal dynamics estimation problem Is, 
given the ''structural*’ Indices of (3-l)-(3-4), i.e. the prespecified basis 


funcclonti f^(x)> eha nimb*r of slgnif leant nodea N, the order of the nodal 
dynamlca n» the nunber of actuators C* and the number of eenaora to 
apply an actual Input aequence at each actuator [u(xj,k)}, neaaure the 
reaultlng output at each aenaor (y(x .k)), and recuraively eatlnate the 

H 

olgenahape parameters c^^ and the modal dynamic parameters a^^ and 
(l.e. NCPf 2n) parameters) to mlnlmlKe the prediction errors 


e(x^,k) - y(Xq,k) - y(Xq,k), (3-5) 

These predicted outputs can be formed from (3-4) with the parameters c^^ 
replaced by their current estimates c^g(h) in 

P 

$.(x,k) - ^ c. (k)f (x) (3-6) 

8-1 “ 

providing the estimated mode shapes used in 

A C ^ 

Uj(k) - I ^.(x. ,k)u(x.,k) (3-7) 

1 j.l 1 3 J 

to provide estimated modal Inputs u^(k) to 

A ® A A A ^ 

yi(h) - bj^4(h)u^(k-l)J (3-8) 

to provide modal output estimates to 

** A ^ 

y(x ,k) » I <).(x ,k)y.(k). (3-9) 

q 1.1 1 q 1 

The approach of [3] and [4], as noted In the preceding section. Is to 
assume that the at® known and from measurements of y(Xq,k-l) to solve 
(3-1) for the yi(k-i) to be used in 


(3-10) 


yj_(V) - I U|^j(k)y^(k-A) + 6^^(k)u^(k-i)J 

1^1 


for comyarlBon with each y^(k) from soluclon of (3-1) given y(x^,k). In order 
for (3-1) to be uniquely solvable the number of sensors Q must equal the 
number of modes N to solve 


‘ y(x2^»k) “ 





y(x2»k) 

• 

■1 

♦^(Xj) *•* 

• • 9 


yjO-) 

• 

9 

y(xQ,k) 


♦ • • 

9 9 9 


9 

y„<« 


(3-11) 


by matrix Inversion or Its equivalent. The question is whether or not such 
a teehnlque can (or should) be incorporated to alter (3-6)- (3-9) and avoid 
the propagation of y^ In (3-8). One problem is the need for as many sensors 
as modes. If this is not feaslblot it seems that the laft-slde of (3-11) 
could be augmented with further measureraants y(x^,k+il). However* this would 
require inclusion of the plant dynamics on the right hand side resulting In 
essentially a multivariable modal observer configuration* which does not 
address the same problem as (3-11) but reverts to the full problem of 
(3-6)- (3-9). (See the next section for consideration of the multivariable 
adaptive observer problem* without the modal structure.) Assuming that 
Q « H Is feasible the problem remains of how to simultaneously provide a 
correction term to (3-6) if the are used In (3-11). That Is* since 


(3-11) with 4*j^(*) replaced by ♦^(•.k) could be used to provide the y^(k-l) 
In (3-8), would y-y using (3-9) provide useful information regarding the 


error in the c^^? 

Since Q =• N Is itself unattractive, consider the apprcach of (3-6)- 
(3-9). First, combine (3-2) and (3-3) to form 




Substituting (3-12) Into (3-1) yields 


^ n Q 

y<x,,k> - ,i /iS)' (k-i) + I <>^(xj )u(xj ,k-i) J). 

j"i (3—1'' 


Then using (3-4) in (3-13) yields 


N P 


y(x .k) - X t 2 c. £ (x )H I e^j^y^(k-i)] 
** 1-1 B-1 ■ ■ ** t-l ^ 


N P 


n CP 


**■ ^ ^ls^s<’‘q>J^ I blit I U <3-14 

1-1 s-1 ■ ‘1 1.1 i*j,i „-i in» ra J 3 


y(x^,k) 


N P n 

iL -L J / li ^^ ia ^ s ^^ q ^ yi ^^"^^ 


1-1 8-1 1-1 


H P n C P 

I I I Z 1 *^ is ^ iro ^ H^s * 

1-1 8-1 1-1 J-1 m-1 IS iro It s q m j j 


(3-15 


Note that even If the measutrement point x^ Is constant over all k, since 

N n n N N 

^ ^ aj^,yi^(k-l) ^ I I I «j'.(x )y.(k-l) 

1-1 ^ ** 1-1 ^ 1-1 i»i i-1 ^ 


n N 

1-1 1-1 ^ 


(3-16 


(3-15) Is not directly transferable to a multi- input (u(Xj,*) over j), 
single-output (y(x^,*) for one q) ARMA process. This Is due to spatial 
coupling, l.e. the output y(x^|^, v) at sensor location x^ Is dependent on the 
past outputs y(x^,») over all q not just those for q = 1. Therefore even 
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th« probictn of estlaaclng the dlfferont pammeter products in 

(3-15) (rather than the desired N(P*f2n) parameters cannot 

be phrased as an equation error parameter estimation [13] probiemy because 
the are not available. (Note for s» m) . Due to the lack of a 

single-output ARMA form an output error formulation [7] [14] is also im- 
possible. One temptation is to relate the fotnt of (3-*S) to that for output 
error estimation. The sieaBurement of y and not y^ requires (3-9) for the 
prediction error and causeB a combination of the y, as In a parallel 
filter implementation, as noted in [15] . The difficulty of (3-6)-(3-9) 
in addition to this structural peculiarity noted in [15] is that the 
effective parameters in this y^ combiner are also now unknown, which leads 
to products of unknown parameters, as is apparent from (3-15). This pro- 
duct form will.be termed (as in [16]) a bllinear-in-the-parametera 


estlmatioa problem, for which there is no known globally stable recursive 
solution. 

Do not confuse the difficulty in establishing a multi-input, single- 
output ARMA model from (3-15) with an. inability to do so in general. 
Clearly if a linear multi-input, multi-output ARMA form exists of the form 


then 


or 


X(k) • I [A.X(k-i) + B.U(k-l)] 
1-1 ^ 


[z“l - I A.z”"^] ^[ I B.z"“^]U(z) 
i-1 ^ i»l ^ 


(3-17) 


(3-18) 


X(z){det(z”l - I A.z""^)}-{Adj[z*^I- f A,z"“^]}{ f B.z""^]U(z) . 
i-1 ^ i«l i-1 


(3-19) 


“9" 


tilling thi Inversi c-traniCorn ylaldt 


XCk) 


S I - 

I a.X(k-i) •¥ I i.U(li-j.), 
1-1 ^ 1-1 ^ 


( 3 - 20 ) 


vheri t - (n) • (dlMnilon ot X). So, a dagrea chango and losa o£ Incarnal 
infomatlon la raqulred to uaa a nodal o£ tba form of (3-20) If It Is usad In 
place of (3*-’17). lliat a couplad nodal, such as (3-17), exists will be tha 
premise for the next two sections. 

Returning to the computation of tha prediction error In (3-5) further 
emphasizes tha blllneac-ln-tha-parametats form of the underlying parameter 
estimation problem. Using (3-13) and the similar form arising from (3-6)- 
(3-9), (3-5) becomes 


e(x^,k) 


I, j l“l. - 


1-1 S «1 


N P 




(3-21) 


Using (3-2) and (3-8) converts (3-21) to 


N P 


)(x ,k) - \ I (c. -c. (k)Jf (x )y.(k) 
1-1 8-1 8 q 1 


1-1 8-1 ‘‘ 1-1 

+{b^^-bj_^(k)}Gi(k-i)] 




t-l 8-1 


i-1 


+ b^j^{u^(k-t)-uj^(k-e)». 


(3-22) 




from (3-3), (3-4), (3-6), and (3-7), (3-22) becomaa 






1-1 4-1 


a-1 


H n 




N P 


n e 


+ I Z I Z I (k)f (x )}f (xju(x.,k-t)j 

1-1 m-1 ^ t-l j-1 ^ 8-1 8 «l a 3 J 


N n 


+ Z Z [yi(k-«)-yi(k-?)][a|^j{ Z 

1-1 1-1 * * ^ g «2 *“ « i 


Due to (3-16) the last term In (3-23) is not n regrossion of o(x , •) that 
could be moved to the left side of (3-23) as in on output error formulation 
[14). Also is unknown In the next to the last term (as In In the 

first term) In (3-23) where It is needed to form the "Input" to the 
^Im^* ‘^dm '’^lm^ segment of the weighted parameter error combination. This 
lattar problem can be solved by approximating b^^^ with (and with y^^) 
as is done in [16] and [17] for a different blllnear-in-the-parameters 
estimation problem. Since for s-m the two terms in c present a 

nontradltlonal problem.) This clearly limits any subsequent estimation 
scheme, based on this approximation, to local convergence. The structure 
of the last term in (3-23) is the more bothersome issue. Assuming that 
^1 ~ ^1 unacceptable method of Ignoring this last term. To assume 

that a^j^ is constant over i is also absurd (unless a breakdown similar 


to (3-'23> la achievable for the plant form In 13«20), which praeently 
seena possible only by losing the decoupled structure of (3-2)}. Another 
improboble situation Is that one set of constant prediction error 


smoothing coefficients would make each of the 1 forward dynamics In an 

r 

output error Identifier error system strictly positive real [18J. The 
form of (3«23) Is enticingly close to a standard output error formulation 
but the problem noted in (3-16) alleviated only by effective solution of 
(3-11) hlndors further consideration of this distributed parameter system 
Identification technique* 

The Justification for developing a atmuJtaneous elgenshape and 
amplitude dynamics estimator is apparent from fl], [4], and the preceding 
section, i.e. the recursive estimation of the b^j, Cj^, and per- 

mits real-time solution of decoupled, Mcalar pole placement problems. This 
contrasts with the large computational effort involved In solution of the 
pole placement problem for a coupled matrix ARMA description as would 
result from an arbitrary fixed choice for the "modal" shape "basis" c^^^. 
Coupling would result in (3-2), l.e. each modal output would be depend- 
ent on past values of all modal outputs and modal inputs, not just Its own 
as In (3-2). Note that "pulse" forms for the c, f (x) products, for 
exi’.uple, would result in a measured, Input (u(Xj, •) - output (y(x^,*)) 
matrix ARMA description, thereby effectively bypassing the modal co- 
ordinate transformations. As this approacli makes the c<!ntrol problem 
solution more involved, the parameter estimation problem becomes solvable. 
Sucli an approach is taken In the next two sect ions. 


4. Simu3.t:9neou« MultivarlabJle Ideatifiaatioii amd Coacyol 

As an alternative to the modal decomiiiOsltiQH repreaentatlon for DPS, 
consider a coupled multi-input multi-output system description r In the 
time domain, this becomes a state space representation with multiple Inputs 
and outputs and suitable state variables. The state space dimension will 
depend on the characteristic behavior of the system In the space and time 
dimensions and the modeling accuracy required. Since DPS require an 
Infinite number of modes In their modal description, this corresponds to 
the need for an infinite dimensional state space for a correspondingly 
complete description. In practice, however, only a finite set of modes 
and therefore finite state space will be assumed for analysis. (See the 
appendix (sec. b) for a more detailed discussion of this modeling issue.) 
Even though a DPS is accurately modeled by a finite state space description, 
this dimension may be too large to raanlpulate in any reasonable real-time 
control application. The required additional reduction in system state 
dimension results in the reduced order control problem, and subsequent ill 
effects caused by modeling Inaccuracy spillover, etc. 112], 

This discussion though, will be limited to the use of simultaneous 
Idontlflcation and control (Indirect adaptive control) on multivariable 
systems without considering these reduced-order modeling effects. The time 
domain approach of simultaneous parameter Identification and state obser- 
vation for use in state variable feedback, will be based on the parameter- 
ized scalar adaptive observer developed by G. Krelselmeier [5] [19]. The 
alternate approach to the Indirect adaptive control problem will be based 
on the frequency domain representation of the multivariable system. This 
is commonly expressed by either a transfer function matrix or matrix 
fraction description relating plant inputs to plant outputs. Here, the 


-13* 


discussion will focus on the frequency domain vet si on of a Loenberger 
obaerver and full atate feedback to accomplish the desired closed loop 
control objective. 

4.1. Time Doealn Approach 

Multivariable extenalon of Kreisselmeier's adaptive obaerver [19] can 

be approached by flrat looking at the major steps In the development for 

single- input single-output systems. If the plant is known to be observable 

and of state dimension Oy l.e. 

X ■ Ax + bu A; nxn, b; nxl, x; nxl 

„ (4-1) 

y « c*x c: nxl 

where (A.c) is observable, an n dimensional observer can be constructed to 
assymptotically estimate the plant states {20] 

X » Fx + g y + h u F; nxn, g; nxl, h: nxl* (4-2) 

For the state estimate error x *» x-x to approach zero, it is required that 
F * A - gc'*’ and b ■ h (4-3) 

where F has eigenvalues strictly in the left half complex plane. This 
observer has the structure in Fig. 4,1. The solution for the state 
estimator is 

X » (g y + h u)dr . (4-4) 

u 0 

However, If the parameters of the plant are unknown i.e. A, B, C are 
known only In dimension, then the observer parameters g and h must be 
estimated such tl-at the state estimate x does indeed converge to the plant 
state X. Kreisselmeier's estimation methods {19] rely on being able to 
separate the observer dynamics and the observer parameters g and h. For 
scalar input and output systems, this is easily accomplished by simply 
commuting terms in the Integrand in (4-4) 

















J y -f h u)dT 

0 


m 


■/*.(A-gcHt-T)y J 1 g ^ rj‘,(A-gcHl:-T) 

0 J LO 


d T 


h 


(4-5) 


The system can now be represented by the structure In PJg, 4-2. 

If this approach la now '» ■■'»? for the tnultl-lnput multi-output case, 
the plant state description bvv < ^ - s 

X w Ax + bu A;nxn li! nxm, u; mxl, xs nxl 

(4-6) 

y » Ox Cs pxn, y; pxl 

where there are m inputs and p outputs. The system is still assumed 
observable. The observer Is described by 

X * Fx -f Gy + Hu P; nxn G: nxp, 11: nxm (4-7) 

with F having strictly left half plane eigenvalues and 

F » A - GC , B * It (4-8) 


to cause assymptotlc state observation. This system has the same structure 
as the scalar system shown in fig. 4-1, except that y and u are now vectors 
and g and h are now matrices. The solution for the state estimate x 
(corresponding to (4-4)) is 

X - + / e^^’'^'^'^^‘^“‘^^(Gy 4- llu)dT (4-9) 

Notice that the outputs y and inputs u will not, In general, commute with 
G and H, respectively. 

To separate the observer dynamics from the obsei'ver parameters con- 
tained In G and II, a unique solution to 

Gy ■ y*G* and Hu ■ u*H* (4-10) 

OR 

^{A-GC)t g . ^(A-GC)*t 

and H - H* (4-11) 

must exist for some set of (■)* quantities. Notice that a unique solution 
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to equation (4-10) will not exist unless G and y span equivalently 
dimensioned suhspaces, i.e. G and y are square. Invertible matrices of 
the same dimension. The same Is true for H and u. The only possibility is 
when y and u are scalars, which reduces the problem to the SISO case 
discussed earlier. 

Equation (4-11) has similar restrictions. Hero, however, If G and H 
are square. Invertible, and of the same dimension as e^^ GC)t^ then unique 
solutions for and H* will exist. For convenience, let ■ ♦(t), then 

_i 

G* * 4>(t)G'{' (t) G{ nxn, Invertible 

_ (4-12) 

H* “ $(t)H$”^(t) H: nxn, invertible 

The solution for the state estimate becomes 

X ■» 4i(t)x« + G* / '}>(t) ^ydT ‘t* H* / ‘{>(t) ^udT . (4-13) 

0 0 

This system has the structure in Fig. 4.3, which is similar to Fig. 4.2. 

The required restrictions for this result, however, are severe! 

(i) The G* and H* estimates must be assymptoticallv invertible for 

A 

X to converge to x . 

(ii) The number of Inputs and outputs must be the same as the number 
of states. 

(ill) The minimum number of states used to describe the system behavior 

must not be overestimated, or Gi'f and H* will never be assymptotical- 

ly Invertible and the state observer may never converge to the 

true plant states. Also, a non-minlmal state description implies 

that some states are either tinohservable or uncontrollable or both. 

If an observed state is uncontrollable, the feedback law may 

require unbounded control inputs In an effort to effect such a 

state, This may drive the -system out of the region of linear 
operation, and Is clearly to be avoided. 
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The epeciellzaticn of the plant model for flexible spacecraft may leaaon 
the severity of the restriction in point (11) » since it may be possible 
to add sensors and actuators to satisfy this point. Restrictions (1) and 
(111) still remain, though, with their inherent numerical problems. 

4.2. Frequency Domain Approach 

Consider the plant having m inputs, p outputs and a pxm proper 


transfer function matrix T(2). Represent T(c) In a left matrix fraction 
description (MFD), not necessarily minimal (8] (l.e. irreducible, relatively 
left prime [20]) 

T(z) » lC^(z)I.(z) (4-14) 

where P(z) and L(z) are polynomial matrices. The elements of F and L are 
polynomials In z whose coefficients are unknown. The system output y(z) 
and input u(z) are then related by 

y(z) * P ^(z)L(z)u(z) (4-15) 

and 


P(z)y(z) ■ L(z)u(z) (4-16) 
Rewrite P(z) and L(z) as sums of products of constant coefficient matrices, 
Pj^ and and powers of z; 


n 


1“0 




y(z) 


V , < 

1“0 


u(z) 


(4-17) 


where n is the largest power of z in P(z) and q is the largest power of z 
in L(z), where q n due to the properness of T(z). Now if P(z) is row 
proper [8], Pq will be invertible* and the taatrix ARMA difference 
equation for y(k) can be given by 


,-l 


r n 

“ I P<y(k-i) + I L.u(k-n+q-i) 
i»l ^ i«0 


y(k) ~ -Q 

Based on this matrix ARMA, some estimation procedure, e.g. [7], can then 


(4-18) 
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be used to estlmete the ARMA coefficients, thereby providing the plant 
parameter estimates for the left MFD. 

A curious quirk in multivariable systems not found In scalar systems 
Is that an equivalent right matrix fraction description for T(s) 

T(z) « R( 2 )Q"^(z) (4-19) 

does not lead, in general, to an ARMA difference equation for y<k). Notice 
that the dual to equation (4-15) 

y(z) « R(z)q"^(z)u(z) (4-20) 

cannot be separated to a dual form of equation (4-16) since R(z) is not 
square unless p « m and even then not necessarily invertible. If such a 
special case holds, then a result similar to (4-16) is 

det|Q(z)]adJ rR(z)]y(z) •• det [R(z) Jadj fQ(z) ]u(z) (4-21) 

If Q(z) and R(z) contains only common unimodular [21] right factors (Q(z) 
and R(z) relatively right prime (r.r.p.) [8J(2l]) then the highest power 
of z on each side of (4-21) will be £ 2n. Here n Is the order of the plant, 
being the number of shifts in a difference equation needed to describe the 
plant. Since an n*"*^ order system has the order of dotjQCz)] equal to n 
(Q(z), R(2) r.r.p.) or greater tijan ii (Q(z)»l<(5!) not r.r.p.) [8,p.l73j and 
the order of adjfR(z)] in z is one or greater, the minimum order on both 
sides of (4-21) Is n. This minimum order of n occurs when R(z) contains 
only constant elements In which case the system has no transmission zeros 
[8,p.l89]. Thus, a minimal right WD can result directly from a matrix 
ARMA difference equation only in the special case where p ■ m and R(z) is 
unimodular. For any particular plant in an adaptive control structure, 
the system order n must also be known for an estimation procedure to 
eventually converge to a minimal right MFD for the plant. However, if 


p m then a right MFD can never result from a matrix ARMA difference 
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Contrast this to the «Ki8tone«t ol tits loft ^SfP basod on a 
matrix ARMA. A left MFD» go dorlvodi always ae^lstS i altbou£h It will be 
dilnlmal only If the orders n and q are known for the plant t therefore, 
hov the estimates of the plant parameters embodied in the left are 
used in feedback control, as veil as the role minimality plays in the 
control effort calculation will now be discussed. 

Using the frequency domain representation of a Luenberger observer 
[8,p.238] provides full plant "state" information for feedback to provide 
arbitrary pole placement< This "transfer function compensation" scheme 


has the structure found in Fig. 4,4* If K(s) and 11 ( 2 ) can be found 

to satisfy the well-known Besout identity [21] 

K(s)q(z) + H(2)R(z) - 1 (4-22) 

then the partial state v can be recreated by measurements of the plant 
Inputs u and outputs y, l.e. 

K(z)Q(z)v + H(2)r(z)v » V (4-23) 

ihe new plant Input is then 

u * r “ F(z)v (4-24) 

and since 

u » Q(z)v (4-25) 

r * lQ(z) + F(z)]v (4-26) 

and with 

y * R(2 >v (4-27) 

y » R(z)lQ(z) + FCzir^r (4-28) 

in which F(z) is chosen such that 

Q(2) + F(z) - Pj(z) (4-29) 


where Pjj(z) is the desired closed-loop denominator matrix. At each Itera- 
tion of an adaptive control algorithm F(z) could be found, based on the 
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currtnt; esKlmata of the plant matrix Q(z), to satisfy (4«>29). Note 
however) that Q(k) ia the denominator matrix of a right MfO» and the 
estimation of the plant parameters produces a left MFD. Horeover) tlie 
right MPD factors Q(*> end R(s) must be relatively right prime for (A-22) 
to be guaranteed a solution for K(z) and H(z) 121). A required step« then) 
will be to find a minimal right Hl^p from a (not necessarily minimal) left 
MFD at each iteration. That this can always be done (but is computation^’ 
ally involved) will now be shown. 

Following tlie procedure in [8M211 for obtaining a greatest common 
left divisor (g.c.l.d.) for the left MFD P""^(z)L(z), find a unimodular 
left multiplier U(z) for the pair P(z) and I.(s) to reduce the construction 
to lower left triangular form 

{F(z),L(z))U(z) « lMCz),0). (4-30) 

The g.c.l.d, M(z) is lower triangular and has the same rank as that of 
{P(z),lCz)l, If (P(z),L(z)] has full row rank m Cor all z, then N(z) has 
rank m for all z, and is therefore uninwdular. In this case, since a 
g.c.l.d. of the pair is unimodular, that pair is relatively left prime. 

If the pair {P(z),L(z)] has rank m for almost, but not all z, then M(z) 
will have similar rank, and M(z) will bo invertible. In this case the 
pair is not r.l.p. but can be made so by eliminating the common non- 
unimodular factor M(z), Multiplying both sides by m"'^(z) yields 

r^(z)[P(z).h<z)lU(z) « 0] (4-31) 

or equivalently (I^ is an mxm identity matrix) 

[P*(z),L*(z))U(z) - (r . 0] (4-32) 

in 

where P*(z) and L*(z) are r.l.p. factors of a left MFD for T(z). Now a 
r.r.p. MFD for T(z) can be obtained by partitioning U(z) as 


K*(ai) R(ss) 

tP*(z),L*(a5)I 

II* (a) -Q(*) 


<A-33) 


• 0 ] . 


W(z) 


Since U(k) Is unimodulert R(x) and Q(ss) are r.r.p. [22] and 

P*(ss)R(z) - L*C*)Q(5t) • 0 (4-34) 

giving 

R(z)Q‘’^(*) - P*"^(z)t*(z) " T(z) (4-35) 

making R(z)Q**^(z) a mlnintal right MFD for T(z). Now since K*(z) and H*(z) 
arc also r»r.p., there exist a K(z) and H(z) that satisfy 

ll(z)K*(z) + K(z)H*(z) ■ (4-36) 

and if K(x) and H(z) are r.i.p. then 

H(z)R(z) - K(s)Q(z) » 0 (4-37) 

resulting In the construction [21 J 


P*Cz) h*(zf 


“k*(z) R(z)“ 

H(z) K(z) 


H*(z) -Q(z) 


U(z) 




(4-38) 


where U(z) is unltnodular» and therefore invertihle, so the solution for H(z> 
and K(z) can be given by 


ni(z),K(z)] - 


0 



U"^(z) . 


(4-39) 


This H(z) and K(z) are a set of polynomial matrices that satisfy equation 
(4-22) for the re-creation of the partial state v. It should be pointed 
out that the above procedure for finding these matrices, while always 
possible, is almost never a trivial matter. The key difficulty is in 
finding the unimodular matrix U(z)» which must be done at every iteration 



of the adaptive aXgorlthM on each new estimate of the P(s) and t(z) 
matrices. 

After the required H(*) and K(z) matrices have been found, a problem 
still exists in the Implementation of the feedback control: H(z), KCz), 
and F(z) are polynomial matrices so that feedback paths are non-proper , 

t 

hence the system is not realizable in real time applications. To over- 
come this difficulty, introduce a stable, Invertible matrix X(z) into 
(4-22) along with the feedback F(z) 


X(z)F(z)K(z)Q(z) + X(z)F(z)H(z)R(z) » X(z)F(z), 


(4-40) 


For any X(z) and F(z), the above calculated H(z) and K(z) are solutions 
to (4-40), and since the r.l.p. pair P*(z) and L*(z) form a left prime 
basis J22) for Q(z) and R-fz), the fenoral solution for (4-40) is 


[ X(z)F(z)K(z), X(z)F(z)G(z) ) 

- X(z)F(z) fK(z) ,H(z) 1 + W(z) lP*(z) ,lXz) ] 


(4-41) 


where W(z) is any polynomial matrix. In this solution, X(z)F(z) [K(z) ,H(z) ] 
Is the particular solution and W(z) [P^(z) ,L>*((z) ] is the homogeneous 
solution. Now choose W(z) such that the general solution 


I X(z)F(z)K(z), X(z)F(z)H(z) ] has row degree less Chan X(z) so that 


X(z)“^ X(z)F(z)K(z) and X(z)“^ 


ZlFlZ 


are proper matrix fractions. 


The system now takes the form shown in Fig. 4-5. The selection of 
W(z) depends on the particular choice of X(z) and the F(z),li(z) and 
P(z) at each iteration of the adaptive algorithm and is not a trivial 
problem. For example, the procedure in [8] requires the Inversion 
of an "ellmlnant matrix” to solve for the X(z)F(z)K(z) and 
X(z)F(z)H(z) (in their notation K(s) and H(s) respectively, which is 


t)rplcally a 10x10 nacrlx for p-m«3. Note also that the matrix manipu- 
lations Involved in any solution technique are complicated by the poly- 
nomial form of the matrix elements, particularly If a machine-calculated 
solution is desired. 

Minimality of the systm MFD Is Important in two respects. The first 
Is that the estimation procedure, used to provide the coefficients of the 
initial left MFD, may require that a minimal structure be known a priori 
for the plant estimation to converse to some meaningful characterization 
of the plant {71. This required a priori Information Includes knowledge 
of the degree n of the plant and effective foreknowledge of the control- 
lability or observability Indices 18], ns well os other structural 
information [23]-r25], such as the relative ayntem degree n-q which is 
related to high frequency behavior. The second reason for minimality 
occurs in the solution for the feedback dynamics baaed on some estimated 
plant MFD. Here, the general solution does not require the initial left 
MFD estimate to be minimal, since the required mintmol right MFD Is found 
in the course of the solution regardless of the minimality of the left MFD. 
However, the solution process is simplified If a minimal left MFD is 
available. 

The key problems with this frequency domain approach are; 

(1) Depending on the particular plant parameter estimation schemes 
and control effort calculation techniques used, minimal system 
descriptions iiwy be required, (This Is analogous to the result 
obtained in the discussion of time domain indirect adaptive control. 

(li) The necessary calculations (left to right MFD conversion and 

compensator parameter calculation) are excessive for performance 
at each step in a real-time adaptive algorithm. 


ThQ next: section will discuss direct adaptive control, which does 
not require explicit identification of the plant parameters, as a possible 
alternative to this scheme and its associated problems. 
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5. )3lract Adaptive Mtiltivaglable Pole Placsmont 

Consider the input , multi-output plant described by the partial 

state description 

A(q'^)s(k) * ii(k) (5-1) 

B(q"^)z(k) - y(k). (5-2) 

“1 -1 

tdiere the nxn polynomial matrix A(q ) in the time delay operator q is 
invertlole such that (5-1) and (5-2) results in the right matrix fraction 
description (MFD) {21,chpt.61 

y(k) - B(q“^A”^(q“bu(k). (5-3) 

where y and u are appropriately dimensioned output and input vectors. 

Kate that a left MFD results froR, a matrix ARMA model as shown in [7] and 
discussed in the preceding section. It turns out that left MPDs are best 
suited for parameter estimation and state observation but right MFDs are 
assumed for feedback control design. As shown In [8] and (9j, in order 
to achieve pole placement via the control law 

C(q"^)u(k) * r(k) + D(q'“^)y(kV, (5-4) 

where C and D are appropriately dimensioned polynomial matrices, C and 
D must be chosen to satisfy 

F(q"b - C(q'^)A(q"^) - D(q"^)B(q"^) , (5-5) 

where F(q~^) is the desired denominator polynomial matrix. This is sub- 
stantiated by substituting (5-1) and (5-2) into (5-4) for 

C(q*'^)A(q”^)z(k) « r(k) + D(tr^)B(q~^)z(k) (5-6) 

or 

[C(q"^)A(q“^) - D(q”^)B(q"^)] 2 (k) - r(k), (5-7) 

Using (5-5) In (5-7) and assuming F is invertible yields 
z(k) * F'‘^(q”^)r(k) . 


(5-8) 
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Uso of (5-2) In (5-8) fesultB In the right MFD 

y(k) « B(q‘^F"^(q”^)r(k), (5-9) 

which In comparison with (5-3) shows that the poles of (5-3) have been 
shifted but that the transmission zeros of (3-3) are unchanged If (5-3) 
and (5-9) are minimal. A critical question is the a priori structuial 
information required to structure C and D and P such that a solution 
exists to (5-5). 

Following the scalar discrete- time strategy fll] (based on the con- 
tinuous-time strategy in (lOJ) for adaptively parameterizing (5-4) without 
a priori specification of A and B in (5-3), a discrete-time multivariable, 
adaptive pole placer will he proposed, Use (5-5) to operate on z yielding 


F(q“^)z(k) » C(q“^)A(q"^) - DCq~^)B(tf bz(k) (5-10) 

If (5-3) is minimal, according to the Bczout identity (21, p. 379] G(q ) 
and H(q ) exist such that 

G(q“bA(q“^) + H(q"^)B(q"h - I. (5-11) 

Inserting (5-11) Into (5-10) yields 

F(q"^)C(q’^)A(q“^z(k) + P(q“^)ll(q“bB(q"^)z(k) 

- C(q“bA(q-bz(k) - l)(q‘^)B(q*’’'-)z(k) . (5-12) 

Using (5-1) and (5-2) in (5-12) yields 

F(q’*^)G(q"^)u(k) + F(q“^)H(q"^)y(k) 

“ C(q~^)u(k) - D(q~^)y(k) (5-13) 

Assuming that F and G (and H) are interchangeable yields 
G(q’^){P(q-hu(k)} + H(q"b{F(q“by(k)} 

- C(q"^){u(k)} + D(q’'^){y(k)} » 0. (5-14) 

As in (10] and (111, estimating G, U, C, and D results In 

e(k) - G(q‘^,k){F(q'’bu(k)} + llCq"^,k){F(q’by(k)} 

- C(q”^,k)(u(k)} + D(q”^,k){y(k)} 
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■ C(q*^,k){u(k)) - D(q”^,k){y(k)} 

M 

- C(q"^»k){F(q"^)u(k)} - fKq“\kHP(q“^)y(k)), (5-15) 

wherep e.g>p 5 ■■ G - G. The error vector e In (5-15) Is recognizable as an 
equation error formulation (13],(26)p which suggests a recursive solution 
of the form 

0(k+l) « 0(k) + P(k)X(k)e(k) (5-16) 


where 

0(k) * (6(q’^,k) f)(q"^,k) G(q"^*k) H(q"^,k)], 


X(k) 


"uXk) 

y(k) 

K(q“bu(k) 
_ P(q“^)y(k)_^ 


(5-18) 


and P is a suitable chosen step-size matrix. Note that this recursion 
could be perforned line by line with each of the entries in the equation 
error vector e, which permits parallel processing thereby reducing the 
computation time per iteration. The number of terms in each of these 
parallel problems increases linearly with the degree of the system in (5-3) 
thereby requiring an increase in the order of the entries in (5-5). 

As noted in [lOj and [llj the stability problem even for the scalar 
case Is unresolved. If the G and H of (5-11) are known exactly and not 
updated In (5-16) then, at least in the scalar case [lOlfll], stability 
can be assured by the technical device of [27] due the atably invertible 
transfer function from u (and y) to e. Comparing (5-5) and (5-11) reveals 
that foreknowledge of G and H is equivalent to foreknowledge of the solution 
to the decoupling, inverse control problem, whicb need not be internally 
stable. Clearly knowledge of this solution corresponds to knowledge of the 
plant parameters. However this encourages the expectation of local 
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stablllty if Che G and H are approximately correct initially. In order to 
retain the almllarity of growth raCea of the input or output and the 
equation error^ |10] suggests bounds on the G and H. How this is to be 
achieved with limited a priori plant information is uncertain; though 
a priori ranges for A and B may translate into acceptable G and H. 

Possibly as uncertainty in A and B increases the acceptable range for G and 
H narrows to Che solution of (5-11) , 

Peculiar to the multivariable case is the structural Information 
required for C and I) to provide a solution to (5-5), especially if F is 
selected in order to form (5-14) from (5-13). In the scalar case, this 
structural information is limited Co plant order and bulk delay (or 
relative degree). The extra complications in the multivariable case, 
just for inverse or model-following control, require foreknowledge of the 
interactor matrix [23] (24] or the Uermltu form [25]. For this pole place- 
ment case, structural constraints may bu difforent. 

The one possibility of this direct adaptive Implementation of pole 
placement, versus indirect schemes, is the seeming possibility of order 
overspecification in the scalar case [tl). This is not possible in the 
indirect case due to the uncontrollable pole-zero cancellation required in 
the identified model for zero identification error. This uncontrollability 
would result In a request for Infinite controller gains leading to adaptive 
controller instability or requiring further logic for avoidance of this 
difficulty. As described in [11] e in the overspecified scalar version 
of (5-15) can be zero with the disappearance of some poles in the overall 
transfer function. This cancellation is stabJ-e 4.ie to the stability of F 
and therefore does not destabilize the adaptive controller. This possibi- 
lity of overspecification in the scalar case raises the hope of overspecified 
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structural Indices In the nultlvarlable case, which could reduce the 
severity of the restrictions mentioned in the preceding paragraph! 
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6. Appendix; Stpayted Variable ModellnR 
BackBround 

A pleueebillty argument will be made for the dlacrete state space 
representation of a distributed parameter system as an approximation to 
the partial differential equation (P.D.E.) representation, subject to a 
limited number of sensor and actuator point locations on ^he system. The 
argument rests heavily on the validity of the separated variable solution 
technique for the P.D.E. . A solution -'mporeci of a factor dependent only 
on time and a factor dependent only on the spatial variables can be 
obtained, provided that the system possesses at least cylindrical symmetry 
about the t-axia in the space spanned by the spatial coordinates and the 
t-coordlnate. Only the class of systems for which this is the case will 
be considered here. Also, the system is assumed linear, 

P.D.E. Representation and Solution 

A linear order P.D.E. in can be represented in general by the 
equation 


L-0 j-0 k-1 




3xj3t 


i 




(6-1) 


I + J - N 


In this equation u(x,t) is a vector of the out-of-equillbrium deflections of 
the system in the spatial coordinates Indexed by k. ja is a function of the 
spatial position vector _x and time. J[(x,t) is a vector forcing function, 
also a function of the position vector x and time. The ^ j^(x, t) terms 
are the coefficients of the various partial derivatives of u. The solution 
of this equation for ii(x,t) is required to satisfy the P.D.E. and be 
uniquely determined by the boundary and initial conditions on some domain 
n in X and t throughout which the P.D.E. representation is valid. 
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If the assumed separated variable solution 

u(x,t> - X(x)T(t) (6-2) 

Is substituted Into the P.D.E.» an Integrating factor can be found such 
that the equation can be arranged having sums of terms, each dependent only 
on X or on t. Those terms that depend on x alone sum to a constant that is 
the negative of the sum of the t-dependent terms. This separation constant 
then appears in the separate solut^’dns for X(x) and T(t), and will be seen 
to play an Important role in the connection between the spatial and temporal 
system solutions. Since the P.D.E. equation is linear, the solution can be 
expressed as the sum of a part due to the natural reaponse to initial con- 
ditions (homogeneous solution) and a part duo to the system forcing function 
(particular solution). The homogeneous solution jj|j(x,t) will be considered 
first. 

The separate homogeneous equations for X(x) and T(t) take the general 

forms 

I B.(T(t).t) [T(t)] » 0 (6-3) 

j»0 ^ df* 

IK 1 

I I y± “ 0 (6-4) 

1*0 k-1 axjj 

where the separation constant a is burled In the 3^ and coefficients. 

J Is the order of the O.D.E. in time, and I is the order of the P.D.E. In 
space. Under certain conditions, the solutions to the above equations 
can be given by linear combinations of orthogonal eigenfunctions 
i/)(t) j|^ (x) 


(6-5) 



T(C) • I 6 (t) 
n-X " " 





(6-6) 


where the conetants A depend on the Initial conditions, and the constants 
are determined by boundary conditions* The eigenfunctions (modes) are 
indexed by n which, ‘along with the separation constant a, determines 
the frequency of the eigenfunctions. Tlius the temporal mode frequencies 
ci>nt snd spatial mode frequencies are related by a. 

The total solutien is then 


i*4l(i£»t) - X(x)T(t) - I 

n*l 

For any particular system represented on the region fi by the above solution, 

the factors ®nd are dependent on the physical nature of the system, 

and the constants depend on the initial conditions of the system. If the 

former factors are known, and the can be sensed or estimated on some 

n 

manner, then the entire status or state of the system is completely known 
in that the output at any time and position can be predicted. 

In practice, the state of the system must be sensed by some finite 
collection of sensors, each of which has a limited area of interaction with 
the system and has a limited frequency response. Therefore, some spatial 
as well as temporal modes will not be sensed. The deflections, velocities, 
etc,, must then be considered as approximations to the true ones at the 
sensor locations. Represent the approximate deflections as made available 
by physical sensors as a finite sum of the eigenfunctions, known constants, 
and constants to be estimated: 
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( 6 - 8 ) 


where cite lowest frequency (spatial and temporal) mode^ are hot neceeasrlly 
the ones sensed » and therefore the index n no longer refers to con- 
secutive mode frequencies* By sensing the u at various points on the 
structurep It la desired that the unknown coefficients in the sum be 
estimated so that the system deflection at any time and at any point In 
space can be predicted. The next section discusses the estimation problem 
for the case of a single point sensor. 


State Estimation - Single Sensor 

If the sensor Is located at some point ^ on the system* the sensed 
deflection at that point is represented by 


I. .... . — — J I.,. , ,..,. 

constants eigenfunctions 

This equation Is a linear corabinatlon of L solutions of the order O.D.E. 
in time. Each solution to the order O.P.K. can be represented by a 
linear combination of solutions of a coupled system of J first order O.D.E.8. 
These J solutions are represented by J state variables. The total repre- 
sentation for u Is then a linear corabinatlon of L sets of J state 
variables. Therefore the output ju can be considered a linear combination 
of L‘J state variables. This can be represented by the following vector- 
matrix equation: 

v(t) « Av(t) ; v(t ) « V 

— •mm 0 ^ 

%(x^,t) * c\(t) 


(6-10) 
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where cha ^ is a vector of the ttate variables. A Is a paraaiecer matrix 
that contains the Information about the natural response to initial con- 
ditions. and ^ is a vector that determines how the state variables com- 
bine to form the output u . Note that the ju vector does not represent 
a vector of outputs but rather a single output on the system at a 
location specified by a position vector. The initial conditions are 
specified by the vector v . It is Important to notice that the expression 
for u in terms of the eigenfunctions In equation (6-9) has been 

replaced by a similar sum of more elementary olgenlunctlons in equation 
(6-10). These elementary functions are all solutions to a first order 
differential equation in time and they all have the form 


^■p(t) - q^exp[rp(t-t^j)l 


V (t) 


( 6 - 11 ) 


where the stero superscript denotes an elementary eigenfunction, the Index 
p ranges from I to L»J, and the and are constants determined by the 
physical properties of the system. The u is then given by: 




L*J 

■« F c V (t) 

p«l ' ' 


( 6 - 12 ) 


where the c are the elements of the c vector. If the physical properties 
P 

A 

of the system are known by some estimation procedure on the sensor output ju, 
and the state variables are known at some time t * through some state 
observation procedure then the state at any time after t * t^^ can be found 
from equation (6-11) and the sensed deflection at any time after t^ can be 
obtained from equation (6-12). It should be pointed out here that some 
vibration modes may not be represented in the sensed deflection of 
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equatlon (6-9) even though their spettei and temporal frequencies as well 
as nodal anplitudes are within the detectable region ol the sensor* This 
is duo to the possibility that the sensor may be located at a point where 
the deflection of tVie body due to some modes Is always too small to detect 
no matter what the inodes' amplitude at other points may be. In this case* 
the sensor is located at a vibrational nodo of those particular spatial 
modes, Such modes are then unobBervabl© in the sensed deflection given by 
equation (6-9), It is important then that the sensor (or sensors) be 
located such that this observability problem does not affect those modes 
of Interest in the system. 

The discussion so far has centered on the deflection of the system 
due to initial conditions only (i,e. homogeneous response). The forced 
(particular) response Involves the additional consideration of external dis- 
turbance forces and actuator forces applied for control purposes. These 
forces can be included in the system model by realizing that external 
forces add energy to each of the characteristic modes In space and time as 
determined by the system's phvslcal properties. With respect to the eigen- 
function expansion description of the system, the forces on the body as 
function of x and time contribute toward spatial modal forces as functions 
of time as expressed by 

m 

- I 4,(t)|„(x) (6-13) 

n*l 

where f_ Is the collective force on the body and are the time varying 

coefficients of the spatial mode shapes ^jj(x). Thus the force on the body 
is represented by a sum of modal forces. For each mode, the deflection 
resulting from a corresponding modal force depends on the physical nature 
of the system such as modal mass, modal damping, modal stiffness, etc. 


i: 

i: 

I 


Juit as in Che expression for the sensed system deflection given by 
equation (6-9) » all physical actuators are limited as to the spatial and 
temporal frequencies they can excite, so the sum In equation (6-13) Is not 
Infinite but limited to say Q excitable modes: 



The modal force emplltudes effectively applied to the system are the Q 
functions of tlme^(t). If these modal force amplitudes are Introduced 
into the O.D.E. in time, equation (6-3), as a non-homogeneous term on the 
right hand side, the result Is a order non-homogeneous differential 
equation. If the equation is linear, each solution can be represented, 
as before, by a linear combination of 3 first order non-homogeneous O.D.E. 
solutions. Using such a set of J equations and corresponding state 
variables for each of the L solutions of equation (6-3) that can be sensed 
by the sensor, a non-honrageneous state space model of dimension L*J is 
obtained 


v(t) • Av(t) + Bj^(x,t) 

" t^v(t) (6-15) 


A 

where B is a matrix whose elements are spatial operators on £(x,t) with 
respect to the elementary eigenfunctions Instead of the original eigen- 
functions as In equation (6-13). These elements depend on the locations 
of the sensors. The system now takes the form of a multi-input, single 
output state space model. Now the system state depends on Initial con- 
ditions as well as the applied forces. Here, a similar problem exists 
with the location of the actuators. Depondlng on the relative location 
of the sensors and actuators, some modes may be excited that have nodes 


e 
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afc the soneoir iocationiit end are therefore unobeervablo. 41so» aome sensed 
nodes may bo Impossible to control if the actuators are located at those 
mode's seros, or nodes. The result in the latter case Is uncontrolXabillty 
, of those system vibrations. For this discussion, It wUl be assumed that 
there are no problems with observability and concroilabitlty in the system 
model . 

State Estimation - Multiple Sensors 

If more than one sensor is used to detect the system state, an 
equation similar to (6-9) can be written for 04ch sensor. The constants 
premultiplying the modes in the sum for each sensor will depend on the 
location necupied on the system by each stMisor. Since the same eigen- 
functions are common to all such sensor deflection representations, the 
same set of state variables can be used to describe the system state at 
any locution, where a differont linear combination of the state variables 
is used at each different location, The vector of sensor measurements is 
given by 

« Cv(t) (6-16) 

where C is a matrix whose i rows £ ^ reflect tlio particular linear com- 
bination of states at each sensor position. 

These particular state variables, being solutions to a modal system 
representation, result In a block diagonal A matrix In equation (6-15). 

Each block represents the solution for each mode, which is ortho. •.onal to 
any other mode, and hence each block is an Independent dynamic system of 
dimension J. 


f 
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$yMWtary 

In more or leee their order of Ifflportence end developnenti the 
aesumptlons and reeulta of this appendix are 

• The P.D«E. representation is linear. 

• this P.D.E. Is separable into apnre and time components. Separate 
apace and time solutions provide that the «'l«onf unction 

form In space does not depend on time .tnd the form In time does 
not depend on the location In Hp.ice wla're monBuraments are taken. 

• The solutions of the separate horaogeneouB space and time equations 
can be represented as infinite auma of ortliugonal eigenfunctions. 

• The frequencies of the space and time eigenfunctions are related 
by the P.D.E. separation constant. 

• Thu constants in the linear combinations of these eigenfunctions 
depend on boundary conditions (for the spaHnl equation) and on 
initial conditions (for the temporal equation) , 

• Thu infinite linear comb Inai ions of eigenfunctions must he con- 
sidered finite for any real liable measurement or actuation due to 
physical limitations. 

• Sensor and actuator placement Is very important with respect to con- 
tro, liability and observability of system modes cjf vibration, 

• A state space representation of the system for arbitrary comblna- 
i lf?ns of actuators and sensors can be theoreticariy found If the 
O.U.K. in time is linear. 

• The order of the state space description depends on the order of 

the 0.0. R. in time and on the number of modes sensable by the sensor. 

• System forces do not alter the eigenfunction forms* but do effect 
Che modal amplitudes, 

• Knowledge of the states in the state variable representation Is 
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mifficient to characterize the 'state of the entire system 
within the accuracy of the sensor measurements. 

• A reduced order model results when the number of states selected 
in the model are fewer than can specify sensor-measurable modes. 
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I • I^ go4<»c t :ion8 

T1j«? three Approaches to adaptive cvintrol of mulci-tnput» multi- 

output, truncated, linear raodela of flexible .structures dlacussod In the 
preceding section (In subsectlunH 3-5) will form the basis of our ongoing 
effortH. Our current empha.sis ta tn reverso order to the order of their 
presentation, l.e. we Intend to investigate, In the following priority 
(1) Direct adaptive multivariable pol« alncement 
(IJ) Simultaneous coupled multivariable systent Identification and 
control via time or frequency domain approaches 
(111) Simultaneous elgenshape and dynamic modal parameter estimation 
and decoupled modal control. 

In order to interrelate these approaclu'S a modeling study Is planned to 
develui> algorithms for eonver.Hion I'row one model form to another, l.e. 
the diHH)upled modal, canonical fU-ite, and matrix fraci on descriptions 
of flexible structures for valruus sensor acuiator locations. Once a 
promising adaptive control candidate emerges It will be Investigated In 
a reduced order setting appropri.ite to DPS or flexible .spacecraft 
control [12]. 
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